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We summarize recent progress in next-to-leading order QCD calculations made using the 
Monte Carlo MCFM. In particular, we focus on the calculations of pp — » Wbb, Zbb and highlight 
the significant corrections to background estimates for Higgs searches in the channels WH 
and ZH at the Tevatron. We also report on the current progress of, and strategies for, the 
calculation of the process pp — > W/Z + 2 jets. 

1 MCFM Background 

With the advent of Run II at the Tevatron we will be able to perform detailed studies of 
femtobarn-level processes for the first time. The production of final states involving heavy 
quarks, leptons and missing energy is particularly interesting since these are common signatures 
for new physics beyond the Standard Model. For example, a light Higgs in the mass range 
110 GeV< mu < 140 GeV predominantly decays into a bb pair. In order to assess search 
strategies and perform meaningful comparisons with the data, a solid knowledge of the Standard 
Model backgrounds is necessary. In an attempt to fill this need, the Monte Carlo program MCFM 
aims to provide a unified description of many of these femtobarn-level processes at next-to- 
leading order (in the strong coupling, as) accuracy. In many cases, the extension to NLO is 
made feasible by the recent calculations of virtual matrix elements involving a vector boson and 
four partons.0 

The philosophy behind MCFM is somewhat similar to that of a general-purpose showering 
Monte Carlo such as Pythiai (although here we are of course working at a fixed order in as)- 
One specifies a set of basic parameters and cuts, selects a process number from a table of included 
processes and then the Monte Carlo MCFM produces a set of relevant distributions. The main pro- 
cesses currently implemented at next-to- leading order are shown in Table |], with various leptonic 
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Figure 1: The m b i distribution of the W (— > u e e + ) bb background to the Higgs signal, at LO and NLO. 



and hadronic decays of the bosons included as further sub-processes. Version 1.0 of the program 



is available for download at http: //www-theory. fnal.gov/people/campbell/mcfm. html. 



Table 1: The main processes included in MCFM at next-to-leading order. 
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2 Higgs Search Using MCFM 

Studies using lowest-order Monte Carlos and other event generators show that for a Higgs in the 
mass range of 100-130 GeV, the most promising channels for discovery at Run II are associated 
Higgs production, a 

pp — ► W(-> ev)H(^ bb), 
pp — ► Z(-» vu,t£)H{-> bb). 

This mass region is particularly interesting in the light of recent hints from LEP2 suggesting a 
Higgs mass mg = 115 GeV. 

Considering first the WH signal, the main backgrounds are, 



pp — ► Wg*(->bb) 
pp — ► WZjf(^bb) 



pp — > t{^bW + )i{^bW~ 
pp — ► W ± *{t{^bW + )b) 
qg — » q't(^bW + )b. 



The signal, Wg* and WZ backgrounds are calculable at the one-loop level in MCFM and the 
remainder at lowest order. In order to study the largest background - the continuum Wg* - 
we introduce a set of standard cuts from the literature, □ use a double 6-tagging efficiency of 
e b i = 0.45 and employ the MRS98 set of parton distribution functions.il When searching for 
the Higgs signal, we are interested in the cross-section as a function of the bb pair mass, m^. 
This distribution is compared at lowest order and next-to-leading order in Figure |l|, with a 
renormalization and factorization scale fi = 100 GeV. From examining this figure, one can see 
that the effect of the next-to-leading corrections is to increase the distribution by a factor of 
approximately 1.5 throughout, thus changing the shape very little. 

The results of a parton-level study for m# = 110 GeV, with all but the small tt and single-top 
backgrounds calculated at next-to- leading order, are summarized in Figure |2[ It is clear from the 




Figure 2: Parton-level signal and backgrounds assuming mn = 110 GeV. The uppermost curve represents the 
signal+backgrounds, with the (almost identical) lower curve representing just the sum of all the backgrounds. 




Figure 3: Sample diagrams contributing to the process pp — > Zbb. Diagram (c) represents a gauge-invariant set 

that does not appear at all in the Wbb calculation. 

figure that the extraction of the Higgs signal requires detailed knowledge of the normalization 
and the kinematics of the backgrounds. 

It is straightforward to perform a similar search in the ZH channel, i Here we shall concen- 
trate on the Z{— * h>v)bb background, calculated at next-to-leading order with MCFM, and compare 
with the Wbb calculation above. The Feynman diagrams required for this calculation - shown 
in Figure |3| - are very similar to the previous case, with additional contributions from gg initial 
states. These gg initial states are important for the Higgs search since they can readily produce 
a bb pair with a large invariant mass. For a conventional scale of 100 GeV, the m b i distribution 
is shown in Figure ||. The effect of the radiative corrections is to produce a large K- factor ~ 1.8 
in the region of interest, m h i from 100-130 GeV. 




3 W + 2 jets: Work in Progress 



The pp — > + 2 jets process can be viewed as an extension of the Wbb and Zbb calculations 
that have already been discussed in the previous sections. There are of course extra parton 
configurations that we must count, but the basic matrix elements (modulo couplings) are the 
same. Specifically, the Wbb calculation contains all the diagrams relevant for qq — > W + q'q'; the 
Zbb process, together with crossings, provides the gg — > W + qq, qq — > W + <?<7, gq — > TV + gg' 
and similar sub-processes. A further complication is that the contribution from the diagrams 
that include real radiation must incorporate the extra singularities due to more instances of 
soft or collinear gluons and collinear quark pairs. In order to simplify the calculation - and to 
reduce the required computation time - we have decided to employ a colour-decomposition of 
the matrix elements. In this procedure the matrix elements are expressed as an expansion in 
1/N (where ./V is the number of colours, 3), with the hope that the sub-leading terms are small 
compared to the leading term. Performing the colour expansion we obtain (V = W ± , Z), 

1 1 

\M NLO (Vqqgg)\ 2 ~ lx^ + pX^ + ^x^ 

\M NLO (VqqQQ)\ 2 ~ - x x Q 3 + S qQ (l x Qo + x 

where the Q, t and Qj represent squared sub-amplitudes for the 2-quark and 4-quark processes 
respectively. Preliminary investigations suggest that the leading piece, Qq is both a good ap- 
proximation to the full matrix element (including both processes) and considerably faster to 
run. For this reason, the inclusion of the term Qq has been the primary focus of the W + 2 jet 
calculation so far. 



4 Conclusions 

We have found large radiative corrections to the Wbb and Zbb processes. These results can 
significantly change estimates of the backgrounds to the processes pp — > WH and pp — > ZH, 
which will be important search channels at the Tevatron. For the more general calculation of 
W/Z + 2 jet production, work is still ongoing but first results should be available soon. 
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